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Abstract 

Phenomenological results of equal time, point to point spatial correlation 
functions of hadronic currents are used to deduce the structure of the QCD 
vacuum. It is found that a model with only quark condensate is not adequate 
to explain the observations. Inclusion of condensate fluctuations (explicit four 
point structure in the vacuum) leads to excellent overall agreement with the 
phenomenological curves and parameters in various channels. 
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The most interesting question in Quantum Chromodynamics (QCD) is concerned with 
the nature of the vacuum state rjTJ . It is well known that the vacuum is non-trivial and is 
composed of gluon and quark field condensates |@||. Much of the understanding of non- 
perturbative phenomena in QCD is expected to result from a proper knowledge of the ground 
state of QCD. 

In this note, we adopt a phenomenological approach to determine the salient features 
of the QCD vacuum. We consider here spacelike separated correlation functions of hadron 
currents in QCD vacuum More precisely, we use phenomenological results of equal time, 
point to point spatial ground state correlation functions of hadronic currents 0] to guide us 
towards a "true" structure of QCD vacuum. 

As a first step, we employ the explicit construct for QCD vacuum with quark and gluon 
condensates proposed by us ||. In this the trial ansatz for the QCD vacuum was 

\vac >= exp (B F * - B F ) exp (B G ] - S G )|0 > (1) 

where |0 > is the perturbative vacuum and B F ^ and B G ^ are the Bogoliubov operators 
corresponding to creation of quark antiquark pairs and gluon pairs respectively. Such a 
construct gives rise to the equal time propagator || as (with x — \x\), 

1 r C(o?),^(0)] \vac} 
J e il *dk [sm2h(k) - (7 ■ k) cos2h{k)] (2) 



Sap{x) = (vac\- 
1 1 



2 (2tt) 3 

where h(k) describes the quark condensate ||. Clearly h(k) gives the free massless 
quark propagator. In the limit of vanishing of constituent quark masses the model predicts 
absence of chiral symmetry breaking. 

We next evaluate the correlation functions of hadronic currents in the vacuum defined 
by Eqs. It turns out that without any approximation, the condensate structure 

uniquely determines the interacting quark propagator and the mesonic(baryonic) correlators 
are essentially squares(cubes) of the propagator. With a Gaussian ansatz for h(k) 0, namely 
sin2/i(/c) = exp(— R 2 k 2 /2) we find that the behaviour of the quark propagator is similar to 



that obtained by Shuryak and Verbaarschott [|7| . Further, the mesonic correlation functions 
were also qualitatively similar to phenomenological correlation functions in all channels 
except for the pseudoscalar(PS) channel. In this channel we did not get the strong attraction 
at intermediate ranges seen in the data. All these results depend very weakly on the choice 
of the functional form of h(k). 

In view of this outcome, it is obvious that some crucial physics is missing from our model, 
and hence the vacuum structure considered by us has to be supplemented by additional 
effects. In our framework, this means that the quark propagators alone do not determine 
the hadronic correlators. This implies that there ought to be explicit contribution arising 
from irreducible four point structure of the vacuum. Alternatively we can describe the 
irreducible four point vacuum structure as a manifestation of condensate fluctuations. 

In order to proceed further with this idea of condensate fluctuations we suggest that the 
normal ordered operators with respect to the vacuum (Eq. (jl|)) do not annihilate the actual 
ground state of QCD which could be more complicated. In this case, we have a more general 
equation 

T^ Q (f)^/3(f)^ 7 (0)?M0) = Sp y (x)Ss a (-x)+ : ^ Q (f)V , /3(x)^ 7 (0)^ (5 (0) : (3) 

where the : denotes normal ordering with respect to the vacuum of Eq.([|). To include the 
effect of fluctuations we may write 

: ^(5)^(^^(0)^(0) := S / 3 7 (f)E <5a (-f) (4) 

so that < vac\HH\vac >= but < fi|££|fi >^ where \Q > is the "new improved" QCD 
ground state including the condensate fluctuations. 

With such a structure for the ground state of QCD the correlator takes the form 

R(x) = Tr [S(x)T'S(-x)T] + Tr [s(f)r'S(-f)r] (5) 

with the extra term arising from the condensate fluctuations. The structure for the fluctu- 
ating field is arbitrary so far with all possible Dirac matrix structures (i.e. 1, 75, 7 M , 7^75, 
Jixu ; H + v ) 



The experimental data demands that we choose the condensate fluctuation field such 
that it contributes maximally in the PS channel and should not affect the other channels 
very much. Such a condition restricts the choice for the fluctuating field to a structure of 
the type 

E aP (g) = + Z%(x) = n\ ( 7 V)a/3 <Hjk + & b aP 0(f) (6) 

where the first term corresponds to vector fluctuations and the second to scalar. ji\ and /i2 in 
the above equations are dimensional parameters which give the strength of the fluctuations 
and 4> k (x) and (f>(x) are vector and scalar fields such that 

< n\(f) k (x)\n >= 0; < n\(f)(x)\n >= o (7) 

and 

<n\(f> i (x)<^(0)\n>=S ii g v (x); < n\(f)(x)(f)(0)\n >= g s (x) (8) 

This implies an approximation that the ground state of QCD is also a condensate in the 
fluctuating fields. The functions gs,v(x) are at this stage arbitrary. It should be pointed 
out that the vector fluctuations contribute only to the PS and nucleon channels, while 
not contributing to the other physical channels. On the other hand the scalar fluctuations 
contribute to all channels but most to the vector and delta channels. It ought to be noted that 
to be consistent with phenomenology of correlators, the contribution of vector fluctuations 
should be greater than that of scalar fluctuations. 

Since we do not expect fluctuations to be important for x — > and for large x , we want 
g (x) to vanish in these two limits. With these properties in mind we take the function g(x) 
as 

9v(x) = [jiiK^x) - ^3^1(^32;)] ; g s (x) = — ^- [^2-^1(^4^) - H5K1 (fax)] (9) 

so that the small x behaviour of the correlator is the same as expected from the asymptotic 
freedom. gv(x) corresponds to taking the fluctuating field condensate function as difference 
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of two propagators with different masses. Such a structure for the condensate function 
arises naturally e.g. in the <j) A field theory in the Gaussian effective potential calculations || 
where \L\ corresponds to the lagrangian mass parameter and [1$ corresponds to the variational 
parameter associated with the Gaussian ansatz. For gs{x) we had to consider a more general 
form (with extra parameters) because the contributions of the propagators and of the scalar 
fluctuations are comparable, requiring delicate cancellations. 

The contributions of the propagator and the fluctuation fields to mesonic and baryonic 
channels are shown separately in cols. 3 and 4 of Table |I[ The parameters of the quark 
condensate function h(k) and the of fluctuating fields gsy are chosen so that our correlation 
functions are similar to those obtained from phenomenology. We have taken R = 0.60/m 
corresponding to < ipip > = —(304.45 MeV) 3 . This value is similar to the value < ipip >= 
— (307.4 MeV) 3 . considered in Ref. PJlfj, though it is higher than the phenomenological 
value. Vector fluctuations dominate in the PS channel beyond x = 0.5/m. This behaviour 
is reproduced by choosing fii = 164.8 and fi^ = 69. In the vector channel the propagator 
and fluctuation contributions are of the same order. If we choose /x 2 — 245, /i 4 = 400, /x 5 = 
1280.77, fiQ = 518 we get reasonable agreement with phenomenological curves of correlators. 
The parameters //i...e are in units of MeV. The resulting correlators are plotted in Fig. |TJ 

We fit our results (solid curves of Fig. |l|) to phenomenologically motivated forms for 
correlators parametrised in terms of mass, coupling and threshold of the corresponding 
particle We used the Marquardt method for non-linear least square fit. The method 
was very stable and the goodness of fit estimated from \ 2 was reasonable for small and 
intermediate x, which is the region characterising the mass, coupling and threshold of a 
general correlator [|J. Our fitted parameters and those similarly obtained in lattice [|J, 
instanton liquid model for QCD vacuum [p|,[10[] and QCD sum rule [nj] calculations are 



tabulated in Table 0. 

As is evident from Fig. [I] and Table [IT] our model of the vacuum gives results for the 
hadronic correlators that are comparable to those in the instanton model H[TIJ and lattice 
calculations [|J . The contribution of condensate fluctuations is most important and in some 



ways is related to the "hidden contribution" discussed by Shuryak [0,[|.To summarise we 
have quite clearly shown that to be consistent with data, QCD vacuum must not only have 
quark condensates but must of necessity also have condensate fluctuations. 
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TABLES 



TABLE I. Meson currents and correlation functions 



CHANNEL CURRENT 



CORRELATION FUNCTIONS 
Without fluctuations 3 " 



R(x) 

Fluctuation contribution 



(Vector(F y ) and Scalar (F s )) 



Vector u'y^d 



[F{x)? + lw e ~ x2/R2 FV = 

F s = 8ir A x 6 g s (2x) 



Pseudoscalar wyz,d 



[F(x)] 2 + f ^e~ x ' 2 / R2 F v = -48TT 4 x 6 g v (2x) 

F s = 4ir 4 x 6 g s (2x) 



Nucleon e abc u a {x)C-i^u b {x) ^j 5 d c (x) [F(x)f + ^^e~ x l R F(x) F v = -^x 6 g v {2x)F{x) 

F s = 2ir 4 x 6 g s (2x)F(x) 



Delta e abc 



u a (x)C lfl u b (x)] u c (x) [F(x)} 3 + l^ e - x2 / R2 F{x) F v = 

F s = 8TT 4 x 6 g s (2x)F(x) 



F(x) 



1 + \x 2 I{x) 



where /(*) = J °° (cos kx - -fM) 



dk 
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TABLE II. Fitted Parameters 



CHANNEL 


SOURCE 




M (GeV) 


A 


V^(GeV) 


Vector 


Ours 




0.78 ± 0.005 


(0.42 ± 0.041 GeV) 2 


2.07 ± 0.02 




Lattice 




0.72 ± 0.06 


(0.41 ± 0.02 GeV) 2 


1.62 ± 0.23 




Instanton 




0.95 ± 0.10 


(0.39 ± 0.02 GeV) 2 


1.50 ± 0.10 




Phenomenolo^ 


;y 


0.78 


(0.409 ± 0.005 GeV) 2 


1.59 ± 0.02 


Pseudoscalar 


Ours 




0.137 ± 0.0001 


(0.475 ± 0.015 GeV) 2 


2.12 ± 0.083 




Lattice 




0.156 ± 0.01 


(0.44 ± 0.01 GeV) 2 


< 1.0 




Instanton 




0.142 ± 0.014 


(0.51 ± 0.02 GeV) 2 


1.36 ± 0.10 




Phenomenolo£ 


;y 


0.138 


(0.480 GeV) 2 


1.30 ± 0.10 


Nucleoli 


Ours 




0.87 ± 0.005 


(0.286 ± 0.041 GeV) 3 


1.91 ± 0.02 




Lattice 




0.95 ± 0.05 


(0.293 ± 0.015 GeV) 3 


< 1.4 




Instanton 




0.96 ± 0.03 


(0.317 ± 0.004 GeV) 3 


1.92 ± 0.05 




Sum rule 




1.02 ± 0.12 


(0.337 ± 0.0.014 GeV) 3 


1.5 




Phenomenolo^ 


;y 


0.939 


? 


1.44 ± 0.04 


Delta 


Ours 




1.52 ± 0.003 


(0.341 ± 0.041 GeV) 3 


3.10 ± 0.008 




Lattice 




1.43 ± 0.08 


(0.326 ± 0.020 GeV) 3 


3.21 ± 0.34 




Instanton 




1.44 ± 0.07 


(0.321 ± 0.016 GeV) 3 


1.96 ± 0.10 




Sum rule 




1.37 ± 0.12 


(0.337 ± 0.014 GeV) 3 


2.1 




Phenomenoloj 


;y 


1.232 


? 


1.96 ± 0.10 
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FIGURES 
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FIG. 1. The ratio of the hadron correlation functions in QCD vacuum to the correlation func- 

R(x) 

tions for noninteracting massless quarks — , vs. distance x (in fm). Our results are given by 

R (x) 

the solid curves. The empirical results determined by dispersion analysis of experimental data in 
Ref.[3] are shown by long dashed lines. The results from lattice calculations and instanton model 
are denoted by dotted and dot-dashed lines respectively. 
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